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Abstract 

We investigated Dirac Particles' Hawking radiation from event horizon of Kerr black hole in 
terms of the tunneling formalism. Applying WKB approximation to the general covariant Dirac 
equation in Kerr spacetime background, we obtain the tunneling probability for fermions and 
Hawking temperature of Kerr black hole. The result obtained by taking the fermion tunneling into 
account is consistent with the previous literatures. 
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Hawking[l] discovered the thermal radiation of a collapsing black hole using the tech- 
niques of quantum field theory in curved spacetime. Since the Hawking radiation relates 
the theory of general relativity with quantum field theory and statistical thermodynamics, 
it is generally believed that a deeper understanding of Hawking radiation may shed some 
lights on seeking the underlying quantum gravity. Since then, several derivations of Hawk- 
ing radiation have been proposed. The original method presented by Hawking is direct but 
complicated to be generalized to other spacetime backgrounds. In recent years, a semi- 
classical derivation of Hawking radiation as a tunneling process Q has been developed and 
has already attracted a lot of attention. In this method, the imaginary part of the action 
is calculated using the null geodesic equation. Zhang and Zhao extended this method to 
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Reissner-Nordstrom black hole[3j and Kerr- Newman black holejjj. M. Angheben et al [5| 
also proposed a derivation of Hawking radiation by calculating the particles' classical action 
from the Hamilton- Jacobi equation, which is an extension of the complex path analysis of T. 
Padmanabhan et al [G] . All of these approaches to tunneling used the fact that the tunneling 
probability for the classically forbidden trajectory from inside to outside the horizon is given 
by 

T = exp(--ImJ) . (1) 

h 

where / is the classical action of the trajectory. The crucial thing in tunneling formalism 
is to calculate the imaginary part of classical action. The difference between these two 
methods consists in how the classical action is calculated. For a detailed comparison of the 
Hamilton- Jacobi method and the Null-Geodesic method, one can see \J\. 

However, most of authors only considered the scalar particles' radiation. Very recently, 
a new calculation concerning fermions' radiation from the stationary spherical symmetric 
black hole was done by R. Kerner and R. B. Mann in ^]. This method have been gener- 



alized to BTZ black hole by us in [9( and dynamical black hole in 10 ]. In this letter, we 



B 

will generalize the tunneling method presented in |8[ to calculate the Dirac particles' Hawk- 
ing radiation from Kerr black hole. Starting with the general covariant Dirac equation in 
curved background, we calculate the tunneling probability and Hawking temperature by us- 
ing WKB approximation. The result obtained by taking the fermion tunneling into account 
is consistent with the previous literatures. 
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The Kerr black hole solution in Boyer-Linquist coordinates is given by 

,2 A - a 2 sin 2 6» , 2 rt . 2/1 r 2 + a 2 - A , , , 
ds 2 = dt 2 - 2asm 2 9 dtdcf) 

(r 2 + a 2 ) 2 - Aa 2 sin 2 # . 2 2 E j2 2 

+ ^— sm 2 6»rf0 2 + — dr 2 + Ed6» 2 , (2) 

2j Za 



where 

E = r 2 + a 2 cos 2 # , 

A = r 2 — 2Mr + a 2 = (r — r + )(r — r_) , 

r± = M ± VM 2 - a 2 . 

We have assumed the non-extremal condition M > a, so that r + and r_ correspond to the 
outer event horizon and the inner event horizon respectively. The determinant of the metric 
is 

y/—g = E sin9 , 

and the inverse of the metric of (t, 0) parts is 

(r 2 + a 2 ) 2 - Aa 2 sin 2 9 



9 tt 



9" 



EA 
A - a 2 sin 2 6» 

EAsin 2 # ' 
a(r 2 + a 2 - A) 
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Now we calculate the Dirac particles' Hawking radiation from the outer event horizon 
of Kerr black hole via tunneling formulism. For simplicity, we only consider the massless 
spinor field \l/ obeys the general covariant Dirac equation 

-^ 7 a <V^ = 0, (3) 

where V M is the spinor covariant derivative defined by = + |<^!f 7[ a 7&] > an d oo a ^ is the 
spin connection, which can be given in terms of the tetrad e£. The 7 matrices are selected 
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as 
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a 2 
\° 2 ° 

where the matrices a k (k = 1, 2, 3) are the Pauli matrices. According to the line element (2), 
the tetrad fields can be selected to be 

T7T/ ' 



= 




e? = 


(0, yf,o,o), 


A 1 


(0,0,-^,0), 


- 

e 3 — 


(0,0,0,-^) . 



V 



exp [-/(*, r,M)] • 



(4) 



We employ the ansatz for the spin- up spinor field \I/ as following j^, [h]] 

/ A(t,r,0,0)\ 


B(*,r,0, 


Note that we will only analysis the spin-up case since the spin-down case is just analogous. 
In order to apply WKB approximation, we can insert the ansatz for spinor field \1/ into the 
general covariant Dirac equation. Dividing by the exponential term and neglecting the terms 
with h, one can arrive at the following four equations 

= «V) + Bj^d r I = , 



/ 



„t<!> 
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(5) 



Note that although A and B are not constant, their derivatives and the components are 
all of the factor h, so can be neglected to the lowest order in WKB approximation. Because 
we only consider the Dirac field outside the event horizon, the condition A > is always 
satisfied in the above equations. The second and fourth equations indicate that 

' d e I + 1-^=8^1 = . (6) 



From the first and third equations one can see that these two equations have a non-trivial 
solution for A and B if and only if the determinant of the coefficient matrix vanishes. Then 
we can get 

(VlFdj - -j^d.i? - ^{d r lf = . (7) 
v-g L 

Because there are two killing vectors (d/dt)^ and (<9/<90) M in Kerr- Newman spacetime, we 
can separate the variables for I(t, r, 9, 0) as following 

I=-Lut + j(f> + R(r,e)+K , (8) 

where u and j are Dirac particle's energy and angular momentum respectively and K is a 
complex constant. Inserting it to equation (7), one can arrive 

n t<f> A 

(v/ 3 ?^ + -j=^3? ~ ^(drR) 2 = . (9) 

V-g L 

Equation (6) indicates that R(r, 9) is a complex function. Now solving the equation (9) for 
definite 9 yields js, 7] 



g u (o y 

dr r.. FZi , ^ 2 A„ 2 „:„2 fl , a(r 2 + a 2 -A) 



± J ^{u^l (r 2 + a 2 ) 2 - Aa 2 sin 2 # - j- 



V(r 2 + a 2 ) 2 - Aa 2 sin 2 o 

The imaginary part of R + can be calculated using the above equation. Integrating the pole 
at the horizon leads to the result (see Q] for a detailed similar process) 

imR± = ± vr(r 2 + a 2 ) ( ^ _ 
r + — r_ 

where Qh = r i ° a 2 is the angular velocity of event horizon. One can see that this result is 
independent of 9. 
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As discussed in the Hamilton- Jacob i method fill. Il2[ |. one solution corresponds Dirac 
particles moving away from the outer event horizon and the other solution corresponds the 
particles moving toward the outer event horizon. The probabilities of crossing the outer 
horizon each way are respectively given by 

2 2 
P out = exp[--ImJ] = exp[--(Imi?+ + ImK)} , 

P in = exp[--Im/] = exp[-|(Imi?„ + ImK)} . (11) 

To ensure that the probability is normalized, we should note that the probability of any 
incoming wave crossing the outer horizon is unity [3]. So we get ImK = — ImiiL. Since 
IulR + = — Imi?_ this implies that the probability of a Dirac particle tunneling from inside 
to outside the evnet horizon is given by 

4 

T = exp[— —lmR + ] , 

= exp[ 4 ; (r + + a2) (c-j^)] , (12) 

where in the last step we set h to unity. It should be noted that the higher terms about 
uj and j are neglected in our derivation and the expression (12) for tunneling probability 
implies the pure thermal radiation. 



From the tunneling probability, the fermionic spectrum of Hawking radiation o 
particles from Kerr black hole can be deduced following the standard arguments 



' Dirac 
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N(w,j) = -z- ( — t^tt , (13) 

where the surface gravity of event horizon. From the tunneling probability 

and radiant spectrum, Hawking temperature of Kerr black hole can be determined as 



r=^ = . (14) 

2tt 4ttM(M + VM 2 - a 2 ) 
In summary, we have calculated the Dirac particles' Hawking radiation from Kerr black 
hole using the tunneling formalism. Starting with Dirac equation, we obtained the radiation 
spectrum and Hawking temperature of Kerr black hole by u sing the WKB approximation. 



The results coincide with the previous literatures 
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